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We report on the first measurement of the triangular U3, quadrangular Vi, and pentagonal us 
charged particle flow in Pb-Pb collisions at y/s^ = 2.76 TeV measured with the ALICE detector 
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of the initial spatial anisotropy and its fluctuations, which provides strong constraints on its origin. 
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as a Mach cone response to fast partons. We show that this structure can be naturally explained 
from the measured anisotropic flow Fourier coefficients. 

PACS numbers: 25.75.Ld, 25.75.Gz, 05.70. Fh 



The quark-gluon plasma (QGP) is a state of matter 
whose existence at high energy density is predicted by 
Quantum Chromodynamics. The creation of this state 
of matter in the laboratory, and the study of its prop- 
erties are the main goals of the ultra-relativistic nuclear 
collision program. One of the experimental observables 
that is sensitive to the properties of this matter is the 



azimuthal distribution of particles in the plane perpen- 
dicular to the beam direction. When nuclei collide at 
non-zero impact parameter (non-central collisions), the 
geometrical overlap region is anisotropic. This initial spa- 
tial asymmetry is converted via multiple collisions into 
an anisotropic momentum distribution of the produced 
particles [T]. 
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The azimuthal anisotropy is usually characterized by 
the Fourier coefficients [2 [3] : 

Vn = (cos[n(0 - *„)]), (1) 

where (j) is the azimuthal angle of the particle, is the 
angle of the initial state spatial plane of symmetry, and n 
is the order of the harmonic. Because the planes of sym- 
metry ^„ are not known experimentally, the anisotropic 
flow coefficients are estimated from measured correlations 
between the observed particles. The second Fourier co- 
efficient V2 is called elliptic flow and has been studied in 
detail in recent years |3] . Large values of elliptic flow at 
the LHC were recently observed by the ALICE collabo- 
ration [S]. 

In a non-central heavy ion collision the beam axis and 
the impact parameter define the reaction plane ^'rp. 
Assuming a smooth matter distribution in the collid- 
ing nuclei, the plane of symmetry is the reaction plane 
= ^'rp and the odd Fourier coefficients are zero by 
symmetry. However due to fluctuations in the matter 
distribution, including contributions from fluctuations in 
the positions of the participating nucleons in the nuclei, 
the plane of symmetry fluctuates event-by-event around 
the reaction plane. This plane of symmetry is determined 
by the participating nucleons and is therefore called the 
participant plane ^'pp [6] . Event-by-event fluctuations of 
the spatial asymmetry generate additional odd harmonic 
symmetry planes \l/„, which are predicted to give rise to 
the odd harmonics like V3 and V5 [7HT5] . 

The large elliptic flow at the Relativistic Heavy Ion 
Colhder (RHIC) [H [H] and at the LHC [5] pro- 
vides compelling evidence for strongly interacting matter 
which appears to behave like an almost perfect (inviscid) 
fluid [12]. Deviations from this ideal case are controlled 
by the ratio 77/s of shear viscosity to entropy density. 
Because the effect of shear viscosity is to dampen all co- 
efficients, with a larger decrease for higher order coeffi- 
cients [TH[T7], it has been argued that the magnitude and 
transverse momentum dependence of the coefficients W3 
and V5 is a more sensitive measure of rj/s [llj . Therefore 
a measurement of these Fourier coefficients at the LHC 
provides strong constraints on the initial geometry, its 
fluctuations, as well as on the shear viscosity to entropy 
density ratio. 

In this paper we report the first measurement of the 
anisotropic flow coefficients ^3, V4, and V5 of charged par- 
ticles in Pb-Pb collisions at the center of mass energy per 
nucleon pair y's^ = 2.76 TeV, with the ALICE detec- 
tor [T5H5D]. The data were recorded in November 2010 
during the flrst run with heavy ions at the LHC. 

For this analysis the ALICE Inner Tracking System 
(ITS) and the Time Projection Chamber (TPC) were 
used to reconstruct charged particle tracks. The VZERO 
counters and the Silicon Pixel Detector (SPD) were used 
for the trigger. The VZERO counters are two scintilla- 
tor arrays providing both amplitude and timing informa- 



tion, covering the pseudorapidity range 2.8 < 77 < 5.1 
(VZERO-A) and -3.7 < 77 < -1.7 (VZERO-C). The 
SPD is the innermost part of the ITS, consisting of two 
cylindrical layers of hybrid silicon pixel assemblies cov- 
ering the range of jTyj < 2.0 and \r]\ < 1.4 for the inner 
and outer layer, respectively. The minimum-bias interac- 
tion trigger required the following three conditions |21) : 
(i) two pixel chip hits in the outer layer of the silicon 
pixel detectors (ii) a signal in VZERO-A (iii) a signal in 
VZERO-C. Deflection of neutral recoils, which is sensi- 
tive to the directed flow of spectators, is measured with 
two neutron Zero Degree Calorimeters (ZDCs) installed 
on each side, 114 meters from the interaction point. Only 
events with a vertex found within 7 cm from the center 
of the detector along the beam line were used in the anal- 
ysis. This is to ensure a uniform acceptance in the cen- 
tral pseudorapidity region \r]\ < 0.8. An event sample of 
5 X 10^ Pb-Pb collisions passed the selection criteria and 
was analyzed as a function of collision centrality, deter- 
mined by cuts on the VZERO multiplicity as described 
in 21 . Based on the strong correlation between the col- 
lision centrality determined by the ZDC, TPC, SPD and 
VZERO detectors, the resolution in centrality is found to 
be < 0.5% RMS for the most central collisions (0-5%), 
increasing towards 2% RMS for peripheral collisions (e.g. 
70-80%). This resolution is also in agreement with our 
Monte Carlo Glauber [32] studies. 

The analysis was, as in [5], performed using tracks mea- 
sured with only the TPC and for tracks using the ITS 
and TPC. These two measurements have very different 
acceptance and efficiency corrections, and provide an es- 
timate of a possible residual bias which may be present, 
after correction, for small values of the harmonics |23j . 
For both measurements, charged particles were selected 
with high reconstruction efficiency and minimal contam- 
ination from photon conversions and secondary charged 
particles produced in the detector material as described 
in ISj. From Monte Carlo simulations of HIJING [21] 
events using a GEANT3 [5S] detector simulation and 
event reconstruction the estimated contamination is less 
than 6% at pt = 0.2 GeV/c and drops below 1% at pt > 1 
GeV/c. In this Letter we present the results obtained 
using the TPC standalone tracks, because of the smaller 
corrections for the azimuthal acceptance. 

We report the anisotropic flow coefficients u„ obtained 
from two-particle correlations and from a four-particle 
cumulant method [26], denoted u„{2} and I'm {4}, respec- 
tively. To calculate the four-particle cumulants we used 
the method proposed in [33]. The f„{2} and w„{4} mea- 
surements have different sensitivity to flow fluctuations 
and contributions from nonflow. The nonflow contribu- 
tion arises from correlations between the particles unre- 
lated to the initial geometry. The contribution from flow 
fluctuations is positive for w„{2} while it is negative for 
'yn{4} [IZ]- Because the odd harmonics are expected to 
be completely due to event-by-event fluctuations in the 
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initial spatial geometry, the comparison of these two and 
four-particle cumulants provides a strong constraint on 
the initial spatial geometry fluctuations. 

The nonflow contribution to the two-particle correla- 
tions is not known and might be significant. We utilize 
four methods to study and correct for nonflow contribu- 
tions to the w„{2} coefficients. First we compare f„{2} 
for like and unlike charge-sign combinations since they 
have different contributions from resonance decay and 
jet fragmentation. Second we used different pseudora- 
pidity gap requirements between the two particles since 
larger gaps reduce the nonflow contributions. Third we 
utilize HIJING (a pQCD inspired model which does not 
include flow) to estimate these contributions and, finally 
we estimate the nonflow from the correlations measured 
in proton-proton collisions. All of these methods indi- 
cate that nonflow effects are smaller than 10%. In this 
Letter we use the dependence of the correlations on pseu- 
dorapidity distance between particles as an estimate of 
nonflow. 

Figure [T^ shows f2, and V4 integrated over the pt 
range 0.2 < pt < 5.0 GeV/c as a function of central- 
ity. The W2{2}, W3{2} and U4{2} are shown for parti- 
cles with |Ary| > 1.0 and corrected for the estimated re- 
maining nonflow contribution based on the correlation 
measured in HIJING. The total systematic uncertainty 
is shown as a band and fully includes this residual cor- 
rection. The measured W3 is smaller than V2 and does 
not depend strongly on centrality. The W3 is compatible 
with predictions for Pb-Pb collisions from a hydrody- 
namic model calculation with Glauber initial conditions 
and rj/s = 0.08 and larger than for MC-KLN CGC ini- 
tial conditions with 7//s — 0.16 llj, suggesting a small 
value of r]/s for the matter created in these collisions. 
The f3{4} is about a factor two smaller than the two- 
particle measurement which can, as explained in [28], be 
understood if U3 originates predominantly from event- 
by-event fluctuations of the initial spatial geometry. For 
these event-by-event fluctuations of the spatial geome- 
try, the symmetry plane ^3 is expected to be uncorre- 
lated (or correlated very weakly [29]) with the reaction 
plane ^'rp, and with 4*2. We evaluate the correlations 
between ^'3 and ^'rp using the first-order event plane 
from the ZDC via v^/^^^ = (cos(3(/)i — 35'rp)) and the 
correlation between 5*3 and 4*2 with a five-particle cor- 
relator (cos(30i -I- 302 — 203 — 204 — 205)) = Vy^^. 
In Fig. [ik i'3/\fnp and Vy^^ are shown as a function of 
centrality. These correlations are indeed, within uncer- 
tainties, consistent with zero as expected from a trian- 
gular flow that originates predominantly from event-by- 
event fluctuations of the initial spatial geometry. 

To investigate the role of viscosity further we calculate 
the ratios V2/S2 and v^/ss, where £2 and £3 are the el- 
lipticity and triangularity of the initial spatial geometry. 




Q 
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FIG. 1. (color online) a) V2, V3 and V4 integrated over the pt 
range 0.2 < pt < 5.0 GeV/c as a function of event centrality, 
with the more central (peripheral) collisions shown on the left- 
(right-)hand side, respectively. Full and open squares show 
U3{2} and «3{4}, respectively. In addition we show v^^^^ and 
173/ ^ J, p, which represent the triangular flow measured relative 
to the second order event plane and the reaction plane, re- 
spectively (for the definitions, see text), b) V2{2, \Ari\ > 1} 
and i'3{2, |Ar7| > 1} divided by the corresponding eccentric- 
ity versus centrality percentile for Glauber [22j and MC-KLN 
CGC [30| initial conditions. 



deflned by: 

(r^cosn(0~ ^„)) 

where the brackets denote an average which tradition- 
ally is taken over the position of participating (wounded) 
nucleons in a Glauber model [22|. 

Under the assumption that Vn is proportional to £„, 
w„{2} is proportional to £„{2} [27]. Figure [i}d shows the 
ratios Vn/Sn for eccentricities calculated with a Glauber 
and a MG-KLN CGC [3D] model, denoted by e^{2} 
and £5^^*^ {2}, respectively. We flnd that for a Glauber 
model the magnitude of f3{2}/£3{2} is smaller than 
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f2{2}/e2{2}, which would indicate significant viscous 
corrections. For MC-KLN CGC calculations the ra- 
tios U2{2}/e2{2} and t;3{2}/£3{2} are almost equal for 
the most central collisions, as expected for an almost 
ideal fluid [TT]. In addition, we notice that the ratio 
'y3{2}/£3{2} decreases faster than i;2{2}/£2{2} toward 
more peripheral collisions which is expected due to larger 
viscous corrections to v^. 

The centrality dependence of the triangular flow dif- 
fers significantly from that of elliptic flow. This might 
be due to two reasons: either the centrality dependence 
of the spatial ellipticity and triangularity are different, 
and/or the viscous effects are different. However, in a 
small centrality range, such as 0-5%, viscous effects do 
not change much and there one might be directly sen- 
sitive to the change in the initial spatial geometry. Our 
calculations show that even in this small centrality range, 
the ratio €2/63 changes significantly which allows us to 
investigate further the geometrical origin of elliptical and 
triangular fiow. In Fig. [2]w2{2} and W3{2} are plotted in 
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FIG. 2. (color online) V2 and 113 as a function of centrality for 
the 5% most central collisions compared to calculations of the 
spatial eccentricities, eJJ'{2} and £^'^'^{2}. The eccentricities 
have been scaled to match the 2-3% data using ki and k2. 

1% centrality bins for the 5% most central collisions. We 
observe that W3{2} does not change much versus central- 
ity (as would be expected if V3 is dominated by event-by- 
event fluctuations of the initial geometry) while, f2{2} 
increases by about 60%. We compare this dependence 
of f„{2} to the centrality dependence of the eccentrici- 
ties £ri{2} for initial conditions from MC-KLN CGC and 
Monte Carlo Glauber. We observe that the weak depen- 
dence of W3{2} is described by both calculations while the 
relative strong dependence of W2{2} on centrality is only 
described for the MC-KLN CGC initial conditions. 

The harmonics V2{'2}, V3{2}, W4{2} and V5{2} as a 
function of transverse momentum are shown for the 30%- 
40%, 0-5%, and 0-2% centrality classes in Fig.jsj For the 
30%-40% centrality class the results are compared to hy- 
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FIG. 3. (color online) V2, W3, V4, vc, as a function of transverse 
momentum and for three event centralities. The full, open 
symbols are for Arj > 0.2 and A77 > 1.0, respectively, a) 
30%-40% compared to hydrodynamic model calculations b) 
0-5% centrality percentile c) 0-2% centrality percentile. 



drodynamic predictions using Glauber initial conditions 
for different values of 77/s [3T]. We observe that, at low- 
Pt, the different pt-dependence of V2 and is described 
well by these hydrodynamic predictions. However, the 
magnitude of V2(j>t) is better described by ij/s = while 
for V3{pt) rj/s = 0.08 provides a better description. We 
anticipate future comparisons utilizing MC-KLN initial 
conditions. 
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For central collisions 0-5% we observe that at pt ~ 2 
GeV/c vj, becomes equal to V2 and at pt ~ 3 GeV/c 
also reaches the same magnitude as V2 and v^. For more 
central collisions 0-2%, we observe that becomes equal 
to V2 at lower pt and reaches significantly larger values 
than V2 at higher-pf The same is true for v^^ compared 

to V2- 

We compare the structures found with azimuthal cor- 
relations between triggered and associated particles to 
those described by the measured Vn components. The 
two-particle azimuthal correlations are measured by cal- 
culating: 



C(A0) 



mixed 



diVsamc/dA0 



(3) 



where : 



d7Vsamc/dA(/) (dA^,„i^ed/dA0) 



is the number of associated particles as function of A0 
within the same (different) event, and A'samc (-^mixcd) 
the total number of associated particles in dA'samo/dA^ 
(dA^niixcd/dA0). Figure [4] shows the azimuthal correla- 

Centrality 0-1%, |ti| < 0.8 

• |Ati| > 1 
V2,3_4,5{2, |Ati|>1} 
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FIG. 4. (color online) The two-particle azimuthal correla- 
tion, measured in < A(j) < tt and shown symmetrized over 
2tt, between a trigger particle with 2 < pt < 3 GeV/c and 
an associated particle with 1 < pt < 2 GeV/c for the 0-1% 
centrality class. The solid red line shows the sum of the mea- 
sured anisotropic flow Fourier coefficients V2, vs, V4 and «5 
(dashed lines). 



tion observed in very central collisions 0-1%, for trigger 
particles in the range 2 < pt < 3 GeV/c with associated 
particles in 1 < pt < 2 GeV/c for pairs in IA77I > 1. 
We observe a clear doubly-peaked correlation structure 
centered opposite to the trigger particle. This feature 
has been observed at lower energies in broader central- 
ity bins [3H [53] , but only after subtraction of the elliptic 
flow component. This two-peak structure has been in- 
terpreted as an indication for various jet-medium modi- 



fications (i.e. Mach cones) [331 [33] and more recently as 
a manifestation of triangular flow [T0Hl3] . We therefore 
compare the azimuthal correlation shape expected from 
V2i vq, and evaluated at corresponding transverse 
momenta with the measured two-particle azimuthal trig- 
gered correlation and find that the combination of these 
harmonics gives a natural description of the observed cor- 
relation structure on the away-side. 

In summary, we have presented the first measurement 
at the LHC of triangular V3, quadrangular Vi and pen- 
tagonal particle fiow V5. We have shown that the trian- 
gular flow and its fluctuations can be understood from 
the initial spatial anisotropy. The transverse momentum 
dependence of V2 and V3 compared to model calculations 
favors a small value of the shear viscosity to entropy ratio 
rj/s. For the 5% most central collisions we have shown 
that V2 rises strongly with centrality in 1% centrality per- 
centiles. The strong change in V2 and the small change 
in V3 as a function of centrality in these 1% centrality 
percentile classes follow the centrality dependence of the 
corresponding spatial anisotropies. The two-particle az- 
imuthal correlation for the 0-1% centrality class exhibits 
a double peak structure around Aip tt (the "away 
side") without the subtraction of elliptic fiow. We have 
shown that the measured anisotropic flow Fourier coeffi- 
cients give a natural description of this structure. 
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